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Abstract: The VuQ model involves the addition of a vector isosinglet up-type 
quark to the standard model. In this model the full CKM quark mixing matrix is 
4x3. Using present flavor-physics data, we perform a fit to this full CKM matrix, 
looking for signals of new physics (NP). We find that the VuQ model is very strongly 
constrained. There are no hints of NP in the CKM matrix, and any VuQ contribu¬ 
tions to loop-level flavor-changing 6 —>■ s, 6 —)■ d and s —?• d transitions are very small. 
There can be significant enhancements of the branching ratios of the flavor-changing 
decays t —)■ uZ and t —)■ cZ, but these are still below present detection levels. 
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1. Introduction 


The standard model (SM) includes three generations of fermions. In particular, 
there are three down-type quarks {Qem = —1/3: d, s, b) and three up-type quarks 
{Qem = 2/3: u, c, t). All quarks with a given charge mix, so that there is a. W 
coupling between each down-type and up-type quark. These couplings are tabulated 
in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. 

Now, there is no a-priori reason for there to be only three down-type and three 
up-type quarks. Indeed, many models of physics beyond the SM include new, exotic 
quarks. The simplest of these consider a fourth generation of quarks (denoted SM4), 
a vector isosinglet down-type quark h' (denoted VdQ; both b'j^ and 6'^ have weak 
isospin / = 0), or a vector isosinglet up-type quark t' (denoted VuQ; both and 
have weak isospin J = 0). 

There are two distinct ways to look for signals of such new physics (NP). The 
hrst is via direct searches at colliders. To date, no signals of exotic quarks have been 
observed; the limits on the masses of these quarks depend on the assumptions about 
how they decay. Some recent results are (this is not exhaustive) */ > 685 GeV 
(95% C.L.) for the SM4 model [|I|, mj,/ > 450 GeV for the VdQ model 0, and mt' > 
687-782 GeV (95% G.L.) for the VuQ model §. 

Second, one can look for indirect signals of the exotic quarks through their loop- 
level contributions to various processes. In fact, it is possible to simultaneously 
consider all such loop-level effects. This is done as follows. Most of these NP ef¬ 
fects are charged-current interactions, which involve the GKM matrix. In the SM, 
the GKM matrix is 3 x 3 and unitary. As such, it is parametrized by four parame¬ 
ters. However, in all NP models the full mixing matrix is larger than 3 x 3, so its 
parametrization requires additional parameters. The idea is then to perform a £t to 
the full GKM matrix using all the data. A signal of the NP will be the non-unitarity 
of the 3x3 GKM matrix. That is, some of the NP parameters will be found to be 
nonzero. 

At hrst glance, the analysis to search for NP is the same for all three models. 
First, in all cases the parametrization of the full GKM matrix has four SM and hve NP 
parameters. Second, one uses the same havor-physics data to perform a combined 
£t to these parameters. This yields the best-£t values of all the parameters, and 
indicates whether any of the NP parameters can be nonzero. However, the key point 
is that the contributions to the havor-physics observables are mo del-dependent. That 
is, the ehects on the observables vary from model to model, so that the analyses are 
not the same for the three models. The SM4 and VdQ models were examined in 
Refs. IQ and |^, respectively. In the present paper we consider the VuQ model 
Bm, in which the full GKM matrix is 4 x 3. 

For the ht, in addition to the six directly-measured magnitudes of GKM matrix 
elements, we include havor-physics observables that have small hadronic uncertain- 
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ties: (i) ex from CP violation in —)■ tttt, (ii) the branching fractions of K~^ —)■ 7r~^uu 
and Kl ^ (hi) Rb and Ai, from Z —)■ bb, (iv) B^-B^ and B^-B^ mixing, (v) the 

time-dependent indirect CP asymmetries in B^ —)■ J/'ipKs and B^ —)■ (vi) 

the measurement of the CP-violating angle 7 of the unitarity triangle from tree-level 
decays, (vii) the branching ratios of the inclusive decays B —)■ XgCC and S —)■ 7 , 
and of the exclusive decay B —)■ iCp+/i“, (viii) many observables in B ^ K*ij,~^n~, 
(ix) the branching ratios of B^ —)■ B^ and B^ —)■ r+z/,-, (x) the 

like-sign dimuon charge asymmetry (xi) the oblique parameters S and T. The 
£t is carried out for rrit' = 800 GeV and 1200 GeV. 

In the VuQ model, the t'^ can mix with the ul, cl and However, because 
the t'j^ and {ul, clRl} have different values of Rl {hi = 0 for t'^^, Rl = ^ for 
{uLyCLyti}), this mixing will induce tree-level Z-mediated flavor-changing neutral 
currents (FGNG’s) among the SM quarks. In particular, this means that D^-D^ 
mixing occurs at tree level. Thus, in principle there can be constraints from the 
experimental measurement of this mixing. Now, in the SM, this mixing is due to a 
box diagram with internal d, s and b quarks. The b contribution suffers a significant 
GKM suppression of 0(A®), so that D^-D^ mixing is dominated by the contributions 
of the internal d and s quarks. Because these quarks are light, there can be large 
long-distance (LD) contributions to the mixing. At present, there is no definitive 
estimate of these LD effects. As a result, we do not have an accurate prediction 
of the value of D^-D^ mixing within the SM, so that this measurement cannot be 
incorporated into the £t. 

Once the £t has been performed, we can then make predictions for other quanti¬ 
ties that are expected to be affected by the t' quark, while still being consistent with 
the above measurements. We examine the following observables: (i) the branching 
fraction of —)■ (ii) the branching fraction of i? —)■ Xgvv, (iii) D^-D^ mixing 

and the branching fraction of —)■ and (iv) the branching fraction of t —)■ qZ 

{q = M,c). 

The paper is organized as follows. In Sec. |^, we dehne the GKM parametrization 
and discuss the measurements used in the fif- The results of the £t are presented 
in Sec. Given these results, we calculate the possible effects of the VuQ model on 
several other flavor observables in Sec. Sec. ^ summarizes the results. 


2. Constraints on the CKM Matrix 


In the VdQ model the GKM matrix is 3 x 4. It was shown in Ref. that this is the 
upper 3x4 submatrix of the 4x4 SM4 GKM matrix, denoted GKM4. Now, there 
are many parametrizations of GKM4. For the VdQ model, it is best to choose one 
in which the new matrix elements Vub', Vcv and Vw take simple forms. With this in 
mind, the Dighe-Kim parametrization of Refs. |^, ^ was used in Ref. 
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WuQ — 


( 2 . 1 ) 


The logic is similar for the VuQ model. In this model the CKM matrix is 4 x 3: 

^ ^ud ^us ^ub 
^cd ^cs ^cb 

Vtd Vts Vtb 

\ Vfd Vfs Vt'b J 

VvuQ is the left-hand 4x3 submatrix of CKM4. Here it is best to choose a parametriza- 
tion of CKM4 in which the new matrix elements Vt'd, Vt's and Vt'b take simple forms. 
We use the Hou-Soni-Steger parametrization |T^, IT^. Here, 




hus — ) 

Vt>d = , W, = 


Vub = , 

Vt<b = -rX , 


( 2 . 2 ) 


where A is the sine of the Cabibbo angle. There are four SM parameters (A, A, C, 
6ub) and hve NP parameters (P, Q, r, 6t'd, St's)- Of the remaining six CKM matrix 
elements, Vud, Vcd and Vcs retain their SM parametrizations: 


Vud = 1 

but Vtdi Vts and Vtb are modihed: 


A 2 \2 

y, W,= -A, W, = l-y, 


(2.3) 


Vu = A\^ (1 - Ce““-) - PrX^e^-'" + -2lCA^e 


iSq, 


SAS, 


Vts = -AX^ - QrX-^e 


+ HAM ^ 


Vtb = 1 


f 2 \ 2 

-r A 
2 


(2.4) 


In the limit P = Q = r = 0, only the elements present in the 3x3 CKM matrix 
retain nontrivial values, and the above expansion corresponds to the Wolfenstein 
parametrization with C = xj and Sub = iwi~^(rj/p). In this limit, Vtb = 1- 
In the VuQ model, r can be nonzero, leading to a deviation of Vtb from 1. 

For the £t, we consider all observables that can constrain the parameters of the 
CKM matrix. The total is written as a function of these parameters, and their 
best-fit values are those that minimize this function. The total function is 
defined as 


2 _ 2 | 2|2 |2 j _2 j _ 2 _|_ 2|2 

Xtotal XcKM ' W Xk^tt+vu ' ' Xz^bb ' SCMh ' Asin2/3 

,2 | 2,2 |2 |2 |2 

Xsin2/3s X'y Xs—tX^ 1+ I— ~r XB—tXg 7 Xb^K /a— Xb—^K* /a— 


+ Xb+ At- + Xb^ 
2 


—AAt+A* 


+ti- + Xb^tu + XaI^ + Aobl: 


ique 


(2.5) 


In our analysis, the x^ of observable A whose measured value is {A^^p ± ) is 

dehned as 


Xa = 


A-A 


exp 


A err 
■^exp 


( 2 . 6 ) 
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iKdl = 0.97425 ± 0.00022 
|K.| = 0.2252 ± 0.0009 
\Vcd\ = 0.230 ±0.011 
|v;,| = 1.006 ±0.023 
\Vub\ = 0.00382 ± 0.00021 
\Vcb\ = (40.9 ± 1.0) X 10-=^ 


21 


7 = ( 68.0 ± 11 . 0 )° 

\ek\ X 10^ = 2.228 ±0.011 
AMd = (0.507 ± 0.004) ps-i 
AM, = (17.72 ± 0.04) ps-^[0 
= 0.00 ± 0.07[iTI 
= 0.68 ± 0 . 02[|1 


^ = 0.00 ± 0.11 
T = 0.02 ±0.12 




B{B X/+r )iow = (1.60 ± 0.48) X lO'^ 

B{B X,£+r)high = (0.57 ± 0.16) X 
lO^GeV^ X (^)(S ^ 77/i+/i")iow = 18.7 ±3.6 
109 GeV^ X {^){B ^ ir/i+/i-)high = 9.5 ± 1.7| 
B{B+ 7r+/iV”) = (2.60 ± 0.61) x 10“^ 
B{K+ TT+z/P) = (1.7 ± 1.1) X 10-1° 
B{Kl /i+/i") = (0 ± 1.56) X 10-9 
B{Bs p+/i”) = (2.9 ± 0.7) X 10-9 1^, 

B{Bd p+p-) = (3.9 ± 1.6) X 10-^9 [g6|, 0 

B{B ^ X, 7 ) = (3.55 ± 0.26) x 10-^ 

B{B ^ r P) = (1.14 ± 0.22) x 10-^ ^ 

= (-4.96 ± 1.69) X 10-9 
7lb = 0.923 ±0.0200 


25 


21 


29 


Rb = 0.2164 ±0.00070 


Table 1: Experimental values of flavor-physics observables used as constraints. For 
Vub we use the weighted average from the inclusive and exclusive semileptonic decays, 
= (44.1 ± 3.1) X 10-^ and = (32.3 ± 3.1) x 10-^. When not explicitly 
stated, we take the inputs from the Particle Data Group 0|. Wherever there are 
asymmetric experimental errors, they are symmetrized by taking the largest side 
error. Also, wherever there is more than one source of uncertainty, the total error is 
obtained by adding these in quadrature. 


In the following subsections, we discuss the various experimental measurements used 
in the fit, and give their individual contributions to Xtotai- 

The current experimental values for the 68 flavor-physics observables enumerated 
in the introduction are listed in Tables ^ and The theoretical expressions for these 
observables require additional inputs in the form of decay constants, bag parameters, 
QGD corrections and other parameters. These are listed in Table p. 


2.1 Direct measurements of the CKM elements 


The latest values for the direct measurements of the magnitudes of the GKM matrix 
elements can be found in Ref. ||0]. The contribution to from these measure¬ 
ments is given by 


2 / |Ks|- 0-2252 W 0.97425 W / |v;,| - 1.006 n 2 

XCKM Q QQQg 0.00022 y ^ V 0.023 ) 

/ |W^|- 0.230 W 0.00382 W 0.0409 W 

^ V 0.011 y ^ V 0.00021 y ^ V 0.001 ) ■ ^ > 

2.2 CP violation in Kl —>■ sk 


In the VuQ model, the mixing amplitude is modified due to an additional con- 
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= 0.1-2 GeV2 

= 2-4.3 GeV2 

= 4.3-8.68 GeV^ 

(^) = (0.60 ±0.10) X 10-7 
(Fl) = 0.37 ±0.11 
(Fi) = -0.19 ±0.40 
(Fa) = 0.03 ±0.15 
(F^) = 0.00 ±0.52 
(F^) = 0.45 ± 0.24 
(F^) = 0.24 ±0.23 
(fJ = -0.12 ±0.56 

(^) = (0.30 ±0.05) X 10-7 
(Fl) = 0.74 ±0.10 
(Fl) = -0.29 ± 0.65 
(Fa) = 0.50 ± 0.08 
(Pi) = 0.74 ± 0.60 
(F^) = 0.29 ± 0.40 
(F^) = -0.15 ±0.38 
(F^) = -0.3 ±0.60 

(^) = (0.49 ±0.08) X 10-7 
(Fl) = 0.57 ±0.08 
(Fl) = 0.36 ± 0.31 
(Fa) = -0.25 ±0.08 
(Pi) = 1.18 ± 0.32 
(Fg') = -0.19 ± 0.16 
(F^) = 0.04 ±0.16 
(F^) = 0.58 ± 0.38 

= 14.18-16 GeV^ 

= 16-19 GeV2 


(^) = (0.56 ±0.10) X 10-7 
(Fl) = 0.33 ± 0.09 
(Fl) = 0.07 ±0.28 
(Fa) = -0.50 ± 0.03 
(R() = -0.18 ± 0.70 
(F^) = -0.79 ± 0.27 
(F^) = 0.18 ±0.25 
(fJ = -0.40 ± 0.60 

(^) = (0.41 ±0.07) X 10-7 
{Fl) = 0.38 ±0.09 
(Fl) = -0.71 ± 0.36 
(Fa) = -0.32 ± 0.08 
(R() = 0.70 ± 0.52 
(F^) = -0.60 ± 0.21 
(F^) = -0.31 ± 0.39 
(F^) = 0.12 ± 0.54 



Table 2: Experimental values of the observables in i? —)■ K* /r’*' ii~ used as con¬ 
straints. These are taken from Refs. [^, Here the errors have been symmetrized 
by taking the largest side error. Also, wherever there is more than one source of 
uncertainty, the total error is obtained by adding these in quadrature. 


tribution coming from a virtual t' quark in the box diagram. There is a sizeable LD 
contribution to the mass difference AM^ in the K system, for which, at present, 
there is no dehnitive estimate. We therefore do not include AM^ in our analysis. 
However, |ei^|, the parameter describing the mixing-induced CP asymmetry in neu¬ 
tral K decays, and which is proportional to Im(Mj^^), is theoretically clean and is 
a well-measured quantity. The theoretical expression for \eK\ in the presence of a t' 
quark is given in Refs. ^ |[. 

To calculate the contribution of \eK\ to xlota\^ quantity 


12\/27r^(AMj^)exp|eE' 

GIM^P^MkBkK 


Im [rfc(V;,V'^)"S(xJ] 


( 2 . 8 ) 


With the experimental and theoretical inputs given in Tables |l| and |], we hnd 


= (1.69 ± 0.05) X 10-' . 


(2.9) 


The QCD correction rjct appears in the theoretical expression of lei^-l. In order to take 
its error into account, we consider it to be a parameter and have added a contribution 
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Gf = 1.16637 X 10-5 Gev-2 
= 0.23116 
oi{Mz) = ^ 
as{Mz) = 0.1184 
mt{mt) = 163 GeV 
mc{mc) = 1.275 ± 0.025 GeV 
mb{mb) = 4.18 ± 0.03 GeV 
Mw = 80.385 GeV 
Mz = 91.1876 GeV 
Mk = 0.497614 GeV 
Mk* = 0.89594 GeV 
Md = 1.86486 GeV 
Mb, = 5.27917 GeV 

Mb, = 5.36677 GeV 
Mb± = 5.27926 GeV 
rrif, = 0.105 GeV 
= 1.77682 GeV 
tb, = (1.519 ± 0.007) ps 


tb, = (1.497 ±0.026) ps 
tb± = (1.641 ± 0.008) ps 
= 0.5765 H 


r]ct = 0.496 ± 0.047 


fK = 0.1561 ±0.0011 16 
Bk = 0.767 ±0.010 H 


AMx = (0.5292 ± 0.0009) x lO-^ps-^ 
Ke = 0.94 ± 0.02 p, II 
K+ = (5.36 ± 0.026) X 10-1^ |3| 

= (2.009 ± 0.017) X 10-^ 


fbd = (190.5 ± 4.2) MeV |T 
fbs = (227.7 ± 4.5) MeV [0 


= (0.216 ± 0.015) GeVil 

^ = 1.268 ±0.063 H 
B{B XJu) = (10.61 ± 0.17) X 10-2 
mc/mb = 0.29 ± 0.02 


Table 3: Decay constants, bag parameters, QGD corrections and other parameters 
used in our analysis. When not explicitly stated, we take the inputs from the Particle 


Data Group |31 


to Xtotai- hold the other QGD correction rjt hxed to its central value because its 
error is very small. The total contribution to Xtotai from \eK\ is then 


Tle/cl 




1.69 X 10 


-7 


0.05 X 10-7 


77^- 0 . 49 6^2 

0.047 


( 2 . 10 ) 


2.3 Branching fraction of the decay 




In Refs. H, it was shown that the LD contribution to B{K^ vr+i/p) is 
suppressed - it is three orders of magnitude smaller than the short-distance (SD) 
contribution. The SM prediction for this observable is therefore under good control. 
The decay —>■ vr+i/P occurs via loops containing virtual heavy particles, and 

hence is sensitive to the t' quark. The theoretical expression for B{K^ tt+z/P) in 
the presence of a t' quark is given in Refs. ii- 

With the inputs given in Tables |] and H, we estimate 


B{K+ 


TT' lyu) 


K+ 


= 3.17 ±2.05 , 


( 2 . 11 ) 
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where 


K+ = rK+- 


3a^B{K^ —)■ TT^e+z/) 


A' 


0 2-4. ■ (2-12) 

Ztt^ sin 

Here rK+ = 0.901 encapsulates the isospin-breaking corrections in relating the branch¬ 
ing ratio of —)■ tt+z/P to that of the well-measured decay —)■ 

In order to include B{K^ vr+z/p) in the £t, we dehne 


XiC+—>.7r+ vv 


[B{K^ —)■ 7r+z/z/)//s;+] — 3.17\2 


2.05 


(2.13) 


2.4 Branching fraction of the decay Kl ^ 

Unlike —)■ tt+z/P, the decay ^ is not cleanly dominated by the SD 

contribution. However, it is possible to estimate the LD contribution to this decay. 
The absorptive LD contribution is estimated using Kl ^ 77 , while the dispersive 
LD contribution is estimated using chiral perturbation theory along with the exper¬ 
imental inputs on various K decays. Due to uncertainties involved in the extraction 
of the dispersive contribution, one can only obtain a conservative upper limit on the 
SD contribution to B{Kl —)■ /r+/i“), which is < 2.5 x 10“® With all the inputs 
given in Tables |l| and ^ we estimate 


B{Kl ^ p+/i-) 


0 ± 0.778 , 


(2.14) 


where 

_ a^B{K^ -)• p+z/^) t{Kl) .8 

^ r(iL+) 


(2.15) 


In the VuQ model, the theoretical expression for B{Ki /i+/i )/is given by 


B{Kl ^ /i+/i-) 




Y{xt) 



(2.16) 


Here Y[x) is the structure function in the t or t' sector |^, ^ 
corresponding structure function in the charm sector, 
value is Pc = 0.115 ± 0.018 [|i 
dehne 


while Pc is the 
Its NNLO QCD-corrected 
In order to include B{Kl P’'"P~) in the £t, we 


Xk 




B{Kl -)■ p+/i )/Kf,-0 


0.778 


+ 


Pc-0.115 
0.018 


(2.17) 


Thus, the error on Pc has been taken into account by considering it to be a parameter 
and adding a contribution to Xtotai- 

















2.5 Z ^ bb decay 


Here we include constraints from Rb and At,, respectively the vertex correction and 
forward-backward asymmetry in Z —)■ bb. The theoretical expressions for Rf, and 
in the VuQ model are given in Ref. |0. We have 


Xz^bb — 


Rfe -0.216 
0.001 


+ 


A - 0.923 \ 2 


0.020 


(2.18) 


2.6 B^-B^ mixing (q = d, s) 

The theoretical expressions for Mf 2 {q = d, s) in the presence of a t' quark, which 
then lead to AM^ and AM^, are given in Refs. @,0. To calculate 
mixing, we use the quantity 


ryd 

■^mix 


GlM^MsAdfl. 

^ d 


With the inputs given in Table |I], we get 


(2.19) 


^ix,exp = (9.12249 ± 1.26905) x 10“^ , (2.20) 


leading to 






9.12249 X 10-^2 


1.26905 X 10-5 


To take Bg-Bg mixing into account, we dehne 


Mr = 


AMgMB^ 1 

AM.Ms^e 


whose measured value is 


( 2 . 21 ) 


( 2 . 22 ) 


Then 


Mr, exp = 21.3831 ± 2.1321. 


Xmh 


/Mr- 21.3831\2 
V 2.1321 ) 


(2.23) 

(2.24) 


2.7 Indirect CP violation in B^ —>■ J/xj^Ks and B^ —)■ J 

The theoretical expressions for Mf 2 (<? = d, s) in the VuQ model are discussed in 
the previous subsection. In the SM, indirect CP violation in Bj^ —)• J/'ipKs and 
—)■ J/'ip4> probes sin 2/3 and sin2/35, respectively. With NP, we have 




JRKs 


Im(M' 


12J 


Mi, 


S 




Im(Mf 2 ) 


M) 


(2.25) 


12 


The experimentally-measured values of sin 2/3 and sin2/3s are given in Ref. Then 

— 0.68\ 2 /Rj/,/,,* — 0.00^ 2 


Xsin 2/3 


0.02 


Asin 2(3s 


J jilx^ 


0.07 


(2.26) 
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2.8 CKM angle 7 

In the Wolfenstein parametrization, the CKM angle 7 = tan“^(7//9), which is the 
argument of Vub- As this angle is measured in tree-level decays, its value is unchanged 
with the addition of a vector isosinglet up-type quark. Therefore the of 7 is given 
by 

'6ub-Q8 (7r/180)\2 

11 (vr/lSO) J 


xl = 


(2.27) 


2.9 Branching ratio of B —>■ Xs I {I = e, fi) 

The quark-level transition b —)• sl^ l~ can occur only at loop level within the SM, 
so that it can be used to test higher-order corrections to the SM, and to constrain 
various NP models. Within the SM, the effective Hamiltonian for this transition can 
be written as 

«.// = -7#E 0((') Oiipi) . ( 2 , 28 ) 

where the form of the operators Oj and the expressions for calculating the coefficients 
Ci are given in Ref. [^|. In the VuQ model only the values of the Wilson coefficients 


^* 7 , 8 , 9,10 are changed via the virtual exchange of the t' quark. The modihed Wilson 
coefficients in the vector-singlet up-quark model can then be written as ii 


cT\^^b) = c,{^,b) + ,, 




(2.29) 


ts^tb 


where j = 7, 8, 9,10. The new wilson coefficients Cj' can be calculated from the 
expression of Cj by replacing rrit by m*/. 

The inclusive decay mode B ^ l~ has relatively small theoretical errors as 
compared to the exclusive decay modes B —>■ {K, K*) l~. However, the inclusive 
decays are less readily accessible experimentally. The branching ratio of R —)■ X* l~ 
has been measured by the Belle and BaBar Collaborations using the sum-of-exclusive 
technique. The latest Belle measurement uses only 25% of its hnal data set [M. The 


BaBar Collaboration has recently published the measurement of B{B —)■ Xg l~) 
using the full data set, which corresponds to 471 x 10® BB events |Q. This is an 
update of their previous result, which was based on a data sample of 89 x 10® BB 
events p8[ . 

The prediction for the branching ratio is relatively cleaner in the low-g^ (1 GeV^ 
< < 6 GeV^) and high-g^ (14.2 GeV^ < < ml) regions. We consider both 

regions in the £t. The theoretical predictions for B{B —>■ Xgl^B) are computed 


using the program Superlso [^, ^ , in which the higher-order and power corrections 
are implemented following Refs. [^1], while the electromagnetic logarithmically- 
enhanced corrections are taken from Refs, 
implemented following Refs. . 


53 . Bremsstrahlung contributions are 
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The contribution to Xtotai 


Xb-i-XsI+I- 


flow 


1.6 X 10-®\2 


0.49 X 10-6 

^ ^B{B Xs 1+ /-)high - 0.57 X 10-6n 2 


(2.30) 


0.23 X 10-6 

where we have added a theoretical error of 7% to B{B —)■ Xg /“)iow, which includes 
corrections due to the renormalization scale and quark masses, and a theoretical 
error of 30% to B{B Xsl~^ l~)high, which includes the non-perturbative QCD 
corrections. 

2.10 Branching ratio of B —>■ Xgj 

The quantity we use for B —)■ 7 is 

5 _ TT f{mc) K{mc) B{B X.q) 

6 a B{B^X,eu,)^ ^ 

where the ratio of the two branching fractions is taken in order to reduce the large 
uncertainties arising from 6 -quark mass. The theoretical expression for B{B Xg 7 ) 
is given in Refs. [^, ^ . From this, one can deduce the expression for R. The measured 
value of R is 

Rexp = 0.1069 ± 0.0120 , (2.32) 

where we have added an overall correction of 5% due to the non-perturbative terms. 
The contribution to Xtota.i 

'B- 0.1069' 


XB^Xs-y — 

2.11 Branching ratio of B 


0.0120 


(2.33) 


K ij,x ^ 


The theoretical expression for {dB/dq^) {B ^ K ^ ) in the SM is given in Refs. 


56| , and can be adapted straightforwardly to the VuQ model. The predictions for 
the branching ratio are relatively cleaner in the low-g^ (1.1 GeV^ < < 6 GeV^) 

and the high-g^ (15 GeV^ < < 22 GeV^) regions. Here, we consider both regions 

in the £t. We use the recent LHGb measurements of {dB/dq^){B K fi~ 


Our analysis of B —)■ X /i+ /r in the low-g^ region is based on QGD factorization 


(QGDf) [^. The factorizable and non-factorizable corrections of 0{as) are included 
in our numerical analysis following Refs. In the high-g^ region, following 


Ref. 1^ , we use the improved Isgur-Wise relation between the form factors which are 
determined using light-cone QGD sum-rule calculations extrapolated to the high-g^ 
region. The contribution to Xtotai from B —)■ iF is 

(^)(B ^ iF/i+/x-)io„ - 18.7 X 10-% 2 


Xb^Ku+ii- — 


6.67 X 10-9 

fr )high — 9-5 X 10 ® 


3.32 X 10-9 


(2.34) 
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where, following Refs. |^, ^], we have included a theoretical error of 30% in both 
low- and high-g^ bins. This is due mainly to uncertainties in the B ^ K form factors. 


2.12 Constraints from B —>■ K* 

The recent LHCb measurements of new angular observables in R —)• K* jj,'^ exhibit 
small tensions with the SM predictions |^, |^. These tensions can be due to NP, but 
can also be attributed to underestimated hadronic power corrections, or can simply 
be a statistical fluctuation. In our analysis, we include all measured observables 
in R —)■ K* fi~ in the low- and high-g^ regions. The experimental results for 
B —)■ K* /i+ decay are given in Table |, and are taken from Refs. [^, Q . 

The complete angular distribution for the decay B —)■ K* /i+ is described by 
four independent kinematic variables: the lepton-pair invariant mass squared g^, two 
polar angles 6^ and 6k, and the angle between the planes of the dimuon and Kn 
decays, 0. The differential decay distribution of R —)■ K* /i"*" can be written as 


d^r[R ^ R:7r)/x+/i-] 

dq^ d cos 9i d cos 9k dcj) 


9 


J 0) • 


(2.35) 


where the angular-dependent term can be written as 


J(g^, 9i,9k,(I)) = Jis sin^ 9k + Jic cos^ 9k + {J2s sin^ 9k + J2c cos^ 9k) cos 29i 
+ J 3 sin^ 9k sin^ 9i cos 20 -|- J 4 sin 29k sin 29i cos 0 

-|- J5 sin 29k sin 9i cos 0 -|- ( Jgs sin^ 9k + Jqc cos^ 9k) cos 9i (2.36) 

-I- Jj sin 29 k sin 9i sin 0 -|- Js sin 29 k sin 29i sin 0 -f Jg sin^ 9 k sin^ 9i sin 20 . 


\L^R /\L,R 


For massless leptons, the Jj’s depend on the six complex K* spin amplitudes Ajj’ , A\_ 
and Aq'^. For example. 


Jis = -A\A\_\^ + \A\? + \Af\^+\A^?]. 


\R\2-\ 


(2.37) 


For massive leptons, the additional amplitude At has to be introduced. In our anal¬ 
ysis, the muon mass is included. 

The analysis of R —?■ K* /i’*' fi~ in the low-g^ region is based on QCDf and 
its quantum held-theoretical formulation, Soft-Collinear Effective Theory (SCET). In 
the limits of a heavy b quark and an energetic K* meson |^, the form factors 

can be expanded in the small ratios Aqcb/ mb and Aqcd/T^, where E is the energy 
of the K* meson. At leading order in l/rub and the seven a-priori independent 
R —)■ K* form factors reduce to two universal form factors H, g 0, n, g. 
The symmetry-breaking corrections of 0{as), both factorizable and non-factorizable, 
are included in our numerical analysis following Ref. |^. Regarding the AQCB/mb 


corrections to the QCDf amplitudes, we do not have any means to calculate them in 
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general. These power corrections can only be estimated by combining QCDf/SCET 
results with a QCD sum rule approach, see Refs. [|^, . 

The analysis of R —)■ K* in the high-g^ region is based on the heavy-quark 

effective theory framework by Grinstein and Pirjol ||^. It was shown in Refs. ^ 
that an operator product expansion is applicable, which allows one to obtain the 
B —)■ R'*/r■*■/!“ matrix elements in a systematic expansion in and in Aqcd/''^ 6- 
The leading Agcn/’^fe corrections are parametrically suppressed and contribute only 
at the few percent level. The improved Isgur-Wise relations between the form factors 
at leading order in 1/m;, lead to simple expressions for the K* spin amplitudes to 

For the form factors in the high-g^ region, we have 


leading order in 1 /m^ |^, ^ 

used the recent lattice results 


Of course, these theoretical predictions have errors associated with them [^, ^ 
f4|, 1^, 1^. The main sources of uncertainties in the low-g^ region, excluding 

uncertainties due to CKM matrix elements, are (i) the form factors, (ii) the unknown 
l/rrib subleading corrections, (iii) the quark masses, and (iv) the renormalization scale 
/ift. In the high-g^ region, there is an additional subleading correction of 0(1/mb) to 
the improved Isgur-Wise form factor relations. For each B —)■ K* observable 

Oj, the theoretical error is incorporated in the fit by multiplying the theoretical 
result by (1 ± Xj), where Xj is the total theoretical error corresponding to the 
observable and can be easily estimated using Table II of Ref. . 

For B —)■ K* jjX jji~^ we use the observables {dB/dq'^), Pi, P 2 , P!i, P^, Pq, Pg 
and Fl in the low-g^ bins 0.1-2 GeV^, 2.0-4.3 GeV^, 4.3-8.68 GeV^, and the high- 

9 1 • -I . 

q 


^ bins 14.18-16 GeV^ and 16-19 GeV^. The SM theoretical expressions for all 


observables in R —)• K* /i“ in the low and high-g^ regions are given in |^, and are 

straightforwardly adapted to the VuQ model by modifying the values of the Wilson 
coefficients as in Fq. (p.29|) . The theoretical predictions for all the B K* 
observables are computed using the program Superlso |^, Q. For each bin, we 
compute the flavor observables and define the 


Xb^k* 




bins 


fi+ fi~ obs.) 


^exp 

3 


Of 


(^i 


2i 


(2.38) 


2.13 Branching ratio of 


77 + fl 




The quark-level transition b —)■ dfi~^ fi~ gives rise to the inclusive semi-leptonic decay 
B^ —>■ Xdji'^ ii~, to exclusive semi-leptonic decays such as B^ fi~, and also 

to the purely leptonic decay B^ —)■ ii~. However, so far, none of these decays have 
been observed. We only have an upper bound on their branching ratios [^, [f^ . 
Recently, LHGb has observed the B~^ 7t~^ ii~ decay with measured branching 

ratio of (2.3 ± 0.6 ± 0.1) x 10“® |^. This is the first measurement of any decay 
channel induced by 6 —)■ d/i’*'p.“. 
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The effective Hamiltonian for the process b ^ d fi~ and the modified Wilson 
coefficients in the VuQ model can be respectively obtained from Eqs. (|2.28|) and 
( p.29|) by replacing s by d. The theoretical expression for —)■ 7r+/i+/r“) is 

given in Ref. ||^. The contribution to Xtotai 


Xb+ 


(B{B^ ^ TT'*' /i’*' /X 


—/i. 




) - 2.3 X 10-®\2 


0.66 X 10-8 


(2.39) 


where, following Ref. ||^, we have included a theoretical error of 10% in B{B^ —)■ 
7T~^ fi~). This is due to uncertainties in the B~^ —)■ tt’*' form factors []8T| . 

2.14 Branching ratio of Bg —>■ (q = s, d) 

The branching ratio of Bg —)■ p” in the VuQ model is given by 

where is defined in Eq. (p.29|) , and (7} 


7}°*’'^ is given by 


T/ T/* 

^tot,d _ ^ I t d t'h 

C'ln — C'lo + -rr^rT~^ 


'tdVtb 


10 • 


In order to include B{Bg fi ) {q = s,d) in the fit, we define 


^lepq 


1671^B{Bq ^ jj,^ jd ) 


GWMB.mlf^^TB, y/l -4(m2/M2J 
Using the inputs given in Tables and ^ we obtain 

-Bieps.exp = 0.025 ± 0.006 , 5iepd,exp = 0.0048 ± 0.0020 . 

The contribution to Xtotai from B{B^ —)■ /i+p”) and B{B^ —)■ id~) is then 


X, Bq — 


B, 


leps 


0.025\2 


0.006 

r u 


+ 


B 


lepd 


0.0020 


0.0048^2 


2.15 Branching ratio of B 

The branching ratio of R ^ r P is given by 


Stt 


m; 

Ml 


fbd\^ub\^'^B±- 


In order to include B{B tv) in the fit, we define 

87iB{B —)■ r P) 


B 


Btau—nu — 


Gp/amJ/^rsVl - m^JAf 
Using the inputs given in Tables and ^ we obtain 

5Btau-nu,exp = (1.779 ± 0.352) X 10"^ 
The contribution to Xtota.i froni B{B v rP) is then 

- 1.779 X 10-%2 


Xb^t V 


0.352 X 10-5 


(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 


(2.45) 


(2.46) 


(2.47) 


(2.48) 
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2.16 Like-sign dimuon charge asymmetry 

The (CP-violating) like-sign dimuon charge asymmetry in the B system is defined 
as 



- N,r- 

Nr + Nr ’ 


where is the number of events of bb —)■ n'^jjL'^X. It can be written as 


(2.49) 


A ^ A d I 

^SL — CsL^SL + C, 


AS 

SL^SL ) 


(2.50) 


where = Im( r ^2 


(9 = s, d), with = 0.594 ± 0.022 and = 0.406 ± 


0.022. The theoretical expression for in the presence of NP is given in Ref. 

has been measured by the D0 Collaboration. The measured value is (—4.96± 


1.53 ± 0.72) X 10-3 
(-2.44 ±0.42) X 10-^ 


29[|. This deviates by 2.7 ct from the SM prediction of A\ 


SL 


IS 


The quantities a, b and c appear in the theoretical expressions for A'gg |^. In 


computing the contribution to from A^gg, one must include the errors in these 
quantities, as well as those in Cgg and Cgg. To do so, we consider all of these as 
parameters and add a contribution to Xtotai- t)e precise. 




4t-(-4.96x 10-3 ) n2 2 

1.69x10-3 ) +Xc. 


(2.51) 


where 


Xc = 


- 0.594n 2 


0.022 

a - 10.5 n2 


± 


- 0.406 


± 


0.022 
5-0.2n2 


1.8 J V 0.1 

2.17 The oblique parameter S and T 


± 


c- (-53.3)\2 
12 ) 


(2.52) 


The theoretical expressions for the oblique parameters S and T in the VuQ model 
are given in Ref. [M. For these non-decoupling corrections we define 


Xoblique 


5-0.0\2 


0.11 


± 


T - 0.02\2 
0.12 ) 


(2.53) 


3. Results of the fit 

We first perform a fit to obtain the Wolfenstein parameters of the standard CKM 
matrix. We then redo the fit, using the theoretical expressions of the VuQ model for 
the observables. We obtain values for the Wolfenstein parameters, as well as for the 
NP magnitudes P, Q and r and the NP phases bt'd and bt's- The results of both fits 
are presented in Table for mt' = 800 GeV and 1200 GeV. 
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Parameter 

SM 

mt>= 800 GeV 

mt'= 1200 GeV 

A 

0.226 ±0.001 

0.226 ±0.001 

0.226 ±0.001 

A 

0.780 ±0.015 

0.770 ±0.019 

0.769 ±0.019 

C 

0.39 ±0.01 

0.44 ± 0.02 

0.43 ±0.02 

^uh 

1.21 ±0.08 

1.13 ±0.11 

1.15 ±0.09 

P 

- 

0.40 ± 0.26 

0.30 ±0.21 

Q 

~ 

0.04 ±0.06 

0.03 ±0.05 

r 

- 

0.45 ± 0.25 

0.36 ±0.22 

St'd 

- 

0.55 ± 0.45 

0.76 ± 0.42 

St's 

- 

0.52 ±3.26 

0.96 ± 1.21 

xVd.o.f. 

71.15/60 

63.35/59 

63.60/59 


Table 4: The results of the hts to the parameters of the CKM matrix in the SM and 
in the VuQ model. 


Quantity 

SM 

mti= 800 GeV 

mti= 1200 GeV 

\Vud\ 

0.9745 ± 0.0002 

0.9745 ± 0.0002 

0.9745 ± 0.0002 

IK.I 

0.226 ±0.001 

0.226 ±0.001 

0.226 ±0.001 

IKbl 

(3.52 ±0.13) X 10-=^ 

(3.92 ±0.24) X 10-3 

(3.85 ±0.21) X 10-3 


0.226 ±0.001 

0.226 ±0.001 

0.226 ±0.001 

IK.I 

0.9745 ± 0.0002 

0.9745 ± 0.0002 

0.9745 ± 0.0002 

IKbl 

0.040 ±0.001 

0.039 ±0.001 

0.039 ±0.001 

\Vtd\ 

0.0084 ±0.0003 

0.0078 ± 0.0005 

0.0080 ± 0.0004 


0.039 ±0.001 

0.039 ±0.001 

0.039 ±0.001 

\Vtb\ 

1 

0.995 ±0.006 

0.997 ±0.004 

\V,'6\ 

- 

0.005 ±0.003 

0.003 ±0.002 

\Vt's\ 

~ 

0.002 ±0.003 

0.001 ±0.002 

\Vt>b\ 

- 

0.101 ±0.056 

0.082 ±0.049 


Table 5: Magnitudes of the 4x3 CKM matrix elements obtained from the £t. 


From Table it can be seen that the three-generation CKM parameters are not 
much affected by the addition of a vector isosinglet up-type quark t'. The allowed 
parameter space for C and Sub expands a little as the constraints on |146| coming 
from the unitarity of the 3x3 CKM matrix are relaxed by the addition of the t' 
quark. The new real parameters, P, Q and r, are consistent with zero. In addition, 
the vanishing of P and Q implies vanishing Vt'd and Vt/si respectively. In this case, 
the phases of these two elements have no signihcance. 

The magnitudes of the elements of the 4x3 CKM matrix, obtained using the 
£t values of Table are given in Table From this Table, we hnd that \Vth\ > 0.98 
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Quantity 

SM 

mt/= 800 GeV 

mt'= 1200 GeV 

\VMVa\ 

0.0084 ±0.0003 

0.0077 ±0.0006 

0.0079 ± 0.0004 

|V,V,;| 

0.0391 ±0.0008 

0.0387 ±0.0011 

0.0386 ±0.001 

\VuV,l\ 

(0.33 ±0.02) X 10-=^ 

(0.30 ±0.02) X 10-3 

(0.30 ±0.02) X 10-3 


- 

(0.47 ±0.40) X 10-3 

(0.28 ±0.26) X 10-3 


- 

(0.19 ±0.32) X 10-3 

(0.12 ±0.20) X 10-3 

|v«v;j 

- 

(0.09 ±0.15) X 10-^ 

(0.05 ±0.09) X 10-^ 


Table 6: In the VuQ model, combinations of CKM matrix elements that control 
mixing and decay in the Bd, Bg and K sectors. 


at 3(j. Now, the direct measurement of \Vtb\i without assuming unitarity, has been 
performed using the single-top-quark production cross section. At the TeVatron one 
hnds \Vtb\ = 1.03 ± 0.06 |^, ^|, while the LHC hnds \Vtb\ = 1.03 ± 0.05 pH ^ 


We therefore see that, although the present direct measurement of \Vtb\ is consistent 
with the SM, a sizeable deviation from its SM value of 1 is not ruled out due to large 
experimental errors. On the other hand, we see that the constraints from present 
flavor-physics data do not allow such a sizeable deviation. We also hnd that the 
allowed values of all of the NP elements of the CKM matrix are consistent with 
zero. Furthermore, the 3cr upper limits on these are |14'd| < 0.01, |Vt's| < 0.01 and 
\Vtib\ < 0.27, indicating that the mixing of t' quark with the other three quarks is 
constrained to be small. 


The values of the magnitudes of the CKM factors that control mixing and decay 
in the Bd, Bg and K sectors are given in Table p. In the 6 —>■ s sector, the NP 
contribution is proportional to the CKM factor Vt'gV^f^. The corresponding CKM 
factor in the SM is VtgV^l- The £t indicates that Thus, the NP 

contribution in the & —)■ s sector is tightly constrained in the VuQ model - large 
deviations from the SM predictions are not possible. This can be seen, for example, 
from the study of the B —)■ K* fi~ observable in the bin 4.3-8.68 GeV^ (see 
Table The disagreement between the experimental measurement of in this bin 
and its SM prediction is around the Aa level. In the SM £t, the Xp' contribution 
to the total Xmin is 16.73, reflecting the large discrepancy between measurement and 
prediction. In the VuQ £t, we find Xp, = 18.18 for mf = 800 GeV (xp, = 17.36 for 
rrif = 1200 GeV), which shows no improvement over the SM. 

The situation is almost the same in the 6 —)■ d and s —)■ d sectors. It can be seen 
from Table | that both |/1| and \VpdVtr^\/\VtdV*g\ are of C>(10-Q. Thus 
the NP contributions in these sectors from the VuQ model are also expected to be 

small. 
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Observable 

Predictions 

SM 

mti= 800 GeV 

mt'= 1200 GeV 

B{Kl —)■ 7r° z/ P) X 10^^ 

2.48 ± 0.29 

3.24 ±0.74 

3.10 ±0.59 

B{B Xgiyp) X 10^ 

2.16 ±0.23 

1.94 ±0.44 

1.95 ±0.40 

Xd 

Unknown 

< 0.08% at 2cr 

< 0.03% at 2a 

B{D /iV) 

^ 3 X 10-^3 

(4.56 ± 10.01) X 10-13 

(1.47 ±2.98) X 10-13 

Bit uZ) 

~ 10-^7 

(1.34 ±2.19) X 10-^ 

(0.50 ±0.89) X 10-^ 

Bit cZ) 

~ 10-^^ 

(1.03 ±2.69) X 10-^ 

(0.39 ± 1.01) X 10-^ 


Table 7: Predictions for observables in the VnQ model. 

4. Predictions for other flavor-physics observables. 


With the constraints fonnd in the previons section for the NP CKM matrix 
elements, it is interesting to see whether any large deviations from the SM are possible 
in other flavor-physics observables. In this section, we provide predictions for some 
of the observables in the VnQ model. These are snmmarized in Table |^. 


4.1 Branching fraction oi Kl ^ i/ i/ 

In the SM, the decay Kl —)■ is dominated by the short-distance loop diagrams 

with top-qnark exchange, while the contribntions dne to the u and c qnarks may be 
neglected. Thus, the t' quark in the loop may give a signihcant contribution. With 
the addition of the t', the branching fraction of Kl —)■ can be written as 0 0 


B{Kl —)■ 7r°i/z/) = Kl 



K{xt) + 



(4.1) 


with 


= (2.31 ± 0.01) X 10“^° , 


r-K^ t{K+) 

The function X{x) {x = relevant for the t and f pieces, is given by 


(4.2) 


X{x) = r]xXo{x) , 


(4.3) 


where 


Xn(x) = 


2 + X 3x — 6 


1 


+ 


X 


(1 


In a; 


X) 


(4.4) 


Above, Tjx is the NLO QCD correction; its value is estimated to be 0.994 |^. rx+ 
summarizes the isospin-breaking corrections in relating K~^ —)■ to —)■ 7r^e~^u, 

while summarizes the isospin breaking corrections in relating Kl ^ 7r°z/n to 
K^ —)■ 
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B{Kl —)■ is a purely CP-violating quantity, i.e., it va ni shes if CP is con¬ 

served. Thus, it is sensitive to non-standard CP-violating phases. Within the SM, 
the branching ratio of —)■ can be predicted with very small uncertainties. 

It is given by 


B{Kl ttVp) = (2.48 ± 0.29) x 10 


-11 


(4.5) 


The main source of uncertainty in the branching ratio prediction is the imaginary 
part of Vtd- Other theoretical uncertainties are less than 2%. Experimentally, this 
decay has yet to be observed. The present upper bound on its branching ratio is 
2.6 X 10“® at 90% C.L. |Q, which is about three orders of magniture above its SM 
prediction. Given the constraints on the 4x3 CKM matrix, the VuQ calculation 
predicts B{Kl —)• = (3.24 ± 0.74) x 10“^^ for mf = 800 GeV ((3.10 ± 0.59) x 

10“^^ for rrit' = 1200 GeV). At 2a, B{Kl —)• < 4.72 x 10“^^, indicating that a 

large enhancement in the branching ratio is not allowed. 

4.2 The branching fraction of .B —>■ Xgiyfy 

In the SM, the decay B —)• is dominated by the penguin and box diagrams 

involving top-quark exchange, and is theoretically clean. Therefore, we expect that 
any additional contributions due to a t' in the loop will be easily identifiable. The 
branching fraction for B —>■ X^z/P in the presence of a t' quark is given by 


B{B X,z/P) = 


o^rjB^B Xceu) 


27r2 sin'^ 6w\Vch\^ f {rhc)n{mc) 


V:,VtsXo{xt] 


1 + 


XQ{Xt' 


v:,ys Xo{xb) 


(4.6) 

The factor fj ^ 0.83 represents the QGD correction to the matrix element of the b —)■ 
sz/P transition due to virtual and bremsstrahlung contributions, /(rhc) is the phase- 
space factor in B{B —>■ XcCp), and k{7jIc) is the 1-loop QGD correction factor. The 
SM prediction for B{B —)■ X^z/p) is (2.16 ±0.23) x 10“®, while in the VuQ model this 
value changes slightly to (1.94±0.44) x 10“^ for = 800 GeV ((1.95±0.40) x 10“^ for 
m*/ = 1200 GeV). Hence a large enhancement of the branching fraction of B —)■ X^z/P 
is not allowed. 

4.3 D^-D^ mixing 

Within the SM, D^-D^ mixing occurs at loop level and involves the lighter quarks 
d, s and b. This implies a strong Glashow-Iliopoulos-Maiani (GIM) cancellation, 
and hence a small SD contribution. Furthermore, the &-quark contribution is highly 
suppressed, 0(A®), so that the mixing is dominated by the d- and s-quark contri¬ 
butions. There are, therefore, large LD contributions to D^-D^ mixing, and indeed 
they dominate over the SD contributions. The present measurement of the D^-D^ 
mixing parameter is 

AMd 


xd = 


D 


= ( 0.8 ± 0 . 1 )%. 


(4.7) 


19 













This is much larger than the short-distance SM prediction. Still, in order to determine 
if the SM can explain this value of xd, one must have an accurate estimate of the 
LD contribution. Unfortunately, this is not available at present. 

As noted in the introduction, the mixing of the with {uL^CL^ti,} will induce 
tree-level Z-mediated FCNC’s among the SM quarks. Thus, in the VuQ model, D^- 
Z)° mixing occurs at tree level. It may therefore provide a much larger contribution 
than that of the (short-distance) SM. Neglecting the SM contributions, in the VuQ 
model D^-D^ mixing is given by 


Xd = 


GF\UucVfhMDBDr{m,,Mz) 


3\/2r 


(4.8) 


D 


where \Uuc\ = Vu^V*^ is the Z-u-c flavor-changing coupling, and r(mc, Mz) = 0.778 is 
the renormalization-group factor. Using /f = 209.2±3.3 MeV |^, = 1.18±0.07 

0 and fD = 0.4101 ps |^, we hnd that, given the constraints on 144^.4, in the 


VuQ model, xd = (0.016 ± 0.034)% for mf = 800 GeV ((0.005 ± 0.010)% for mf = 
1200 GeV). Thus at 2cr, xd < 0.08%. We therefore see that the SD contribution in 
the VuQ model falls far below the observed value of D^-D^ mixing. 

4.4 Branching fraction of —)■ 

Unlike D^-D^ mixing, the SM prediction for the branching fraction of —)■ fjGiji~ 

can be estimated fairly accurately, even after including the LD contribution. The 
SM prediction for the —)■ branching ratio is ~ 3 x 10“^^, hence highly 

suppressed. Thus, —)■ has the potential for large NP contributions. At 

present, we only have an experimental upper bound on the branching ratio: B{D^ —)■ 
< 7.6 X 10“® at 95% G.L. [Q, which is several orders of magnitude larger 
than the SM prediction. 

Within the VuQ model, —)■ occurs at tree level due to Z-mediated 

FGNG’s. Neglecting the SM contribution, the branching ratio in the VuQ model is 
given by 


/i V) = 


_ GpralflMD 


327rF D 


Ami 

^-\U 
\ ^ \ 


ml 


UC 


(4.9) 


For mp = 800 GeV, B{D^ pV) = (4-56 ± 10.01) x lO-^^ ((1.47 ± 2.98) x lO-^^ 
for mp = 1200 GeV). Thus, at 2(j, B{D^ —)■ < 2.46 x 10“^^. We therefore 

observe that the branching ratio of —)■ can be enhanced by an order of 

magnitude above its SM value, but this is still far below the present detection level. 

4.5 Branching fraction of t —>■ qZ (q = c, u) 

Within the SM, the branching ratios of the FGNG top decays t uZ and t ^ cZ are 
~ 10“^^ and ~ 10“^*^, respectively |^, The present upper bound on B(t qZ) 
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is 0.21% at 95% C.L. |Q. The discovery potential of B{t —)■ qZ) is ~ 10 -10 ° at 


ATLAS and CMS. The SM valne of B(t ^ qZ) is thus far below the detection level 
for these decays. This implies that these decays can only be observed if NP enhances 
their branching ratios by many orders of magnitude above their SM values. 

This may be possible within the VuQ model, as here, due to Z-mediated FCNC’s, 
these decays occur at tree level. Neglecting the SM contribution, the decay rate for 
t —)■ qZ is given by 


qZ) = 
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32 sin^ 6w cos^ 6w ' 

where rrit = 173.2 ± 0.9 GeV E3 and \Uqt \ = VqiV^^. As Vtb in this model is close to 
unity, we can approximate the top width by r(t —)■ hW^), which at leading order is 
given by 


T{t bW^) = 


a 


-\Vtb\^ 


m 


3 r 


M'i- 


16 sin^ 9w ' 

The branching ratio of t —)■ qZ is therefore given by 
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3—+ 2—^ 


B{t qZ) = (0.463 ±0.001) 
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(4.12) 


Using the values of parameters given in Table ||, we obtain \Uut\ = (0.53±0.43) x 10 ^ 
((0.33±0.29)X10-3) and |Urf| = (0.47±0.61)xl0-3 ((0.29±0.37)xlO-^) for mj! =800 
GeV (1200 GeV). This leads to B{t uZ) = (1.34±2.19) x 10“^ ((0.50±0.89) x 10“^) 
and B{t ^ cZ) = (1.03 ± 2.69) x 10"^ ((0.39 ± 1.01) x 10"^) for mt> = 800 GeV 
(1200 GeV). Therefore, the FGNG branching ratios can indeed be enhanced by many 
orders of magnitude above their SM values. However, they are still two orders of 
magnitude below the present detection level for these decays. 


5. Conclusions 

In this paper we consider the VuQ model, in which a vector isosinglet up-type quark 
t' is added to the standard model (SM). In the VuQ model, the full GKM quark 
mixing matrix is 4 x 3, and is parametrized by four SM and hve new-physics (NP) 
parameters. The NP parameters include three magnitudes and two (GP-violating) 
phases. We perform a £t using flavor-physics data to constrain all GKM parameters. 
The purpose is to determine whether there are any indications of NP, such as the 
non-unitarity of the 3x3 SM GKM matrix, or, equivalently, nonzero values for 
some of the NP parameters. And even if there is no evidence of NP, we would 
like to ascertain whether sizeable NP effects are still possible in other flavor-physics 
observables, while being consistent with the constraints found in the £t. 

The £t involves 68 flavor-physics observables. No evidence for NP is found: the 
values of the three NP magnitudes are consistent with zero, in which case the two 
NP phases have no signihcance. Specihc results include the following: 
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The deviations of the CKM matrix elements Vts and Vtd from their SM predic¬ 
tion are small. 


• At 3a, \Vtb\ > 0.98. Any large deviation of \Vtb\ from unity is therefore not 
possible in the VuQ model. 

• The 3a upper limits on the new elements of the VuQ CKM matrix are: \Vtid\ < 
0.01, \Vt's\ < 0.01 and \Vt'b\ < 0.27, indicating that the mixing of t' quark with 
the other three quarks is constrained to be small. 

Turning to possible NP effects in the VuQ model, we hnd that any NP contri¬ 
butions to 6 —)■ s, 6 —)■ d and s —)■ d transitions are tightly constrained. We also 
hnd, 

• A large enhancement of SD contribution to Xd (i.e., D^-D^ mixing) is not 
allowed. 

• The branching ratio of —)■ can be enhanced by an order of magnitude 

above its SM value, but this is still far below the present detection level. 

• The branching ratios of the havor-changing decays t —)■ qZ [q = c, u) can be 
enhanced by many orders of magnitude. However, they are still two orders of 
magnitude below the present detection level. 

In summary, current havor data puts extremely stringent constraints on the VuQ 
model. There are no hints of NP in the CKM matrix. Furthermore, the ht to the 
data indicates that any VuQ contributions to loop-level havor-changing 6 —)■ s, 6 —)■ d 
and s —)■ d transitions are very small. There can be signihcant enhancements of the 
branching ratios of t —)■ uZ and t ^ cZ decays, but these are still below detection 
levels. 
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